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We focus on very general, very large U, doped Mott Insulators with arbitrary hopping and inter- 
actions. We provide simple testimony to the competition between magnetic and superconducting 
orders in these systems. By mapping hard core bosons, spinless, and spinful fermions onto XXZ 
models, we aim to make very simple precise statements. We try to address optimal and expected 
filling fractions of holes within the plane and on stripes in a variety of hole and hole pair geometries. 
We examine the role of attractions/repulsion amongst hole pairs and single holes, and provide trivial 
expected numerical values for filling fractions in various scenarios. We demnostrate that plaquette 
states seem to naturally provide the correct stripe filling fractions. 



I. INTRODUCTION AND OUTLINE 

Throughout this paper we aim to address questions 
concerning planar and linear (stripe) filling fractions of 
doped Mott insulators by trying to see how much may be 
gleaned by mapping such systems onto spin models. The 
answers that we find depend greatly on the the assumed 
form of our constituent particles (single holes, pairs) and 
their various geometries (plaquettes, rungs, bonds). We 
will demonstrate that plaquette states seem to naturally 
provide the correct stripe filling fractions. By mapping 
onto spin models, we will also be able to examine various 
competing orders. Amongst other things, we will demon- 
strate (not unexpectedly) that in the limit of large on- 
site repulsions U, magnetic and superconducting orders 
always compete. 

The outline of the paper is as follows. In Section(II), 
we lay out the general standard model of Doped Mott 
Insulator which will form the focus of our discussion. 
In Section(III), we examine the problem of spinless and 
spinful particles in the plane by a mapping to a spin 
model. We will see, in a special yet general set of mod- 
els another rigorous example of the competition between 
magnetic and superconducting orders. Within these spin 
models that result for large U, we illustrate how mag- 
netic order as well as the physically transparent number 
order (both portrayed by S z in two different spin repre- 
sentation) compete with the dual superconducting order 
encapsulated by bilinears in the planar S± multiplied by 
exponentials in of a topological nature. 

In Section(IV), we map a model of hard core bosons 
(hole pairs) in non-overlapping plaquettes with arbitrary 
finite range interactions and hoppings onto a planar XXZ 
model. We find that if the physics is indeed dominated 
by such attractive pair states then, within the ground 
state, the average hole occupancy per site within the op- 
timal doped state may be 1/4. The large deviation from 
the observed optimal doping in most doped Mott insu- 
lators is hardly surprising and points to the inadequa- 
cies of looking at attractive plaquette pairs alone within 



the plane: not unexpectedly, an analysis of the cuprates 
based solely on notions of doped Mott insulators and con- 
currently assuming only Bose pairs of the plaquette size 
might be flawed. 

We next examine, in Section(V), plaquette pairs on 
bond-centered stripes and find, under a variety of cir- 
cumstances, near 1/4-filling. We discuss the possibility 
of phase separation or modulation of the hole pairs along 
the bond centered stripe. When we relax the condition 
that diagonal pairs do not have to be in non-overlapping 
plaquettes, we will find, in the large hopping amplitude 
(t) limit, a 1/3-filling fraction of stripes. Higher and simi- 
lar filling fractions are found for pairs on single rungs and 
legs. From this simple exercise we conclude that bond 
centered stripes cannot, perhaps, be described by simple 
pairs roaming the stripe axis, and that if pairs indeed do 
dominate the asymptotic low energy stripe scale physics 
as we suggested in an earlier paper [1] , then they must be 
effectively confined to non-overlapping plaquettes in ac- 
cord with pictures suggested by the DMRG calculations 
of White and Scalapino [2] . 

We the proceed with a similar analysis of fermions on 
bond centered stripes (Section(VII)), and find that these 
lend themselves to near 1/4 filling. Our description for- 
tifies earlier work by Nayak and Wilczek [3] . 



II. THE MODEL 

We will focus our attention on the relatively standard 
model of doped Mott insulator, the extended Hubbard 
model: 

h = - *ij( c L c J<r + c L c j>) + Yl Un ^ n ^ 

(ij)o i 

where c ia creates an electron on site i with spin a and 
j is a nearest neighbor of i. This model contains both 
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the movement of the electrons (hopping) (t, kinetic en- 
ergy) and the interactions of the electrons if they are on 
the same site (U, potential energy). The Mott insulat- 
ing nature is captured by this on-site repulsion which 
greatly inhibits hole motion. We have added an addi- 
tional term representing all possible number-number in- 
teractions (of all ranges); these may be result from myr- 
iad interactions- e.g. Coulomb repulsions, interactions 
mediated by phonons. The number occupancy 

T^ia Cecity 

ni = ^n trT . (2) 

(7 

As well known, in the infinite U limit, trivial spin- 
charge separation occurs in any dimension. The charge 
degrees of freedom may be trivially encapsulated by spin- 
less degrees of freedom. For a review of this principle, the 
reader is invited to read Appendix(A). 

At large U, the extended Hubbard model of Eqn.(l) 
may be related to an extended t-J model, 

H = - M C L C ^ + H.c.) + J ijSi ■ Sj 

(ij) 
ij 

where Si = ^2 aa i c\<j® <?■<?' c i-<?' is the spin of the electron at 
site i, a are the Pauli matrices, and there is a constraint 
of no double occupancy of any site i (rii = J^u c l a°i,<y 
has expectation values or 1). The reduction to a 
Hilbert space where no doubly occupied sites occur (the 
Gutzwiller projection) will be automatically incorporated 
in all things to come. 

III. COMPETING ORDERS FOR FERMIONS IN 
THE PLANE 

A. Spinless Fermions- A Competition Between 
Charge And Superconducting Orders 

With possible applications to the limit of infinite on- 
site repulsion (U) in mind, we now examine the compe- 
tition of charge and superconducting order within the 
plane. All that we will detail below can be derived 
straightforwardly. We take a slightly longer route in or- 
der to hightlight the simple similarities between spin and 
charge when looked at through the prism of the Jordan- 
Wigncr representation of S=l and S=l/2 problems re- 
spectively. Although not of any use for most practical 
applications, a high dimensional Jordan Wigner trans- 
formation has been devised by Fradkin [4], and later ex- 
tended by Eliczer and Semcnoff [5]. Very nice novel ex- 
tensions to various spinful cases were recently advanced 
by Batista and Ortiz [6]. The basic message is that the 



classic Jordan- Wigner [7] transformation rigidly linking 
spinless fermions and S=l/2 spins in one dimension can, 
quite naturally, be extended to higher dimensions. The 
only complication is that now the kink operators that 
code for the statistics transmutations become high di- 
mensional topological objects. To be more precise, the 
standard one dimensional string operator appearing in 
the usual Jordan- Wigner transformation is replaced by 
its more general counterpart 

^•=exp[i^0(fc,J)ng], (4) 
k 

with 9{k,j) the angle between (k — j) and a fixed ray (in 
the one dimensional case the "angle of site" 9 reduces to 
either it (for j < k) or zero (for k > j) leading to the 
standard string operator). In terms of these new kink 
operators, the Jordan- Wigner transformation then reads 

St = c] Kj , 

s; = K] Cj , 

S*= nj -l/2. (5) 

The simple pair operator Ajj = ejej may be expressed in 
terms of the spin variables. In the aftermath, we find that 
the superconducting pairing operator may be expressed 
as a product of kink variables with the XY components 
of the spins: Sf 1 and S^. The incompatibility of charge 
(n or S z ) and superconducting phase (descendant from 
S±) orders is trivially reflected from the non- vanishing 
commutator 

[S Z ,S ± ]^±S ± . (6) 

Explicitly, inverting all matters to the spin representa- 
tion, 

4 = S+K- 1 = S+ cxpH E ^ , (?) 

f 

the commutator 

[n k , Aij] = (S ik + 5 jk )Aij, (8) 

which is compatible with Ay = exp[i(0p + <f>f)], and 
n k = —ig^u with C denoting charge. 

If Eqn.(l) captured all of the physics, then computing 
optimal doping would amount to the determination of the 
value of (S z ) when we aim to maximize the topological 
pairing order parameter A. 

B. Spinful Fermions- A Competition between 
Magnetic and Superconducting Orders 

Just as in the spinless case, the simple commutation 
relations that we will soon derive can be immediately 
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seen by direct computation (no Jordan- Wigner transfor- 
mations are neccessary). Nevertheless, in order to high- 
light the similiraties between the charge and spin sectors 
as doublet (S=l/2) and triplet (S=l) representations of 
similar entites we will employ the Jordan- Wigner repre- 
sentation once again. Batista and Ortiz [6] extended the 
Jordan- Wigner transformations to spinful fermions. The 
operator 

c] a = 4(1 - ni- a ), (9) 

along with its conjugate, may be transformed into spin 
S=l operators in much the same way as for the spinlcss 
case. We note that within the large U limit, cj CT — ► c] CT . 
Any Hamiltonian in {ct CT , Cj CT } has the Gutzwiller (no 
double occupancy) projection automatically built into it. 
The general t-J Hamiltonian of Eqn.(3) will undergo no 
change when expressed in terms of {c] CT ,Cj CT } instead of 

The transverse (XY) components of the spin may be 
written as [6] 

S+ = 2 1 ' 2 {z) >a=+ K i + Kfe ia= -). (10) 

The z-component of the spin transforms as the on-site 
magnetization, 

Sj z = Tlj t<7 =+ — Tlj t<7 ——, (11) 

with nja = c^Cja. As before, 
K = cxp[i£ 6(k,j)n k ] = cxp[i^0(fe, 3)(Stf], (12) 

k k 

where the, new, second equality follows from Eqn.(ll) 
and the inbuilt constraint of no double occupancy. When 
we employ the inverse transformation, we find that Ay = 
c ia c j-a contains, similar to the case before, the trans- 
verse components of the spin S± multiplying topological 
operators K~ x , which are simple exponentiated products 
in {(S? ) 2 }. Once again the non-commuting character of 
the spin components (Eqn.(6)) disallows concurrent ideal 
magnetic and superconducting orders. Explicitly, 

[S' k ,A ij ] = ~{5 ik -5 jk )A ij . (13) 
This is compatible with 



Ay=exp[i(^ (14) 
with the superscript M denoting "Magnetic" . Trivially, 
[S z (q),A(p)]^0, (15) 
with q and p arbitrary wavenumbers. 



Examining Eqs.(8,13), we note that, apart from an im- 
portant minus sign, magnetism competes with supercon- 
ductivity in much the same way as the charge (number) 
field competes with pairing for spinlcss fermions. This 
is, of course, not all too surprising. We merely wished 
to highlight the similarities between the spin and charge 
sectors as S=l and S=l/2 representations of a similar 
problem when viewed through the Jordan- Wigner trans- 
formation. 

All of what we detailed henceforth was for singlet pair- 
ing. For a hypothetical triplet pair operator A^ = 

CfoCja the commutation relations are trivially identical to 
those of the number and pairing correlations of Eqn.(8), 

[n k ,^] = \{5 ik + 5 jk )A^. (16) 

To summarize, in a well known simple rubric, 

[Spin, Pairing] <~ [Charge, Pairing] ^ 0. (17) 

We further reiterate that the charge and spin sectors can 
be viewed as a doublet (S=l/2) and tripet (S=l) rep- 
resentations of similar entities. The spin component S z 
takes on different roles the two different Jordan- Wigner 
transformations. This is simply yet another way of view- 
ing matters. The similarity between the two sectors is 
highlighted in the SO(5) theory of S. C. Zhang [8]. 

IV. PLANAR PLAQUETTE PAIR STATES 

Much of our approach henceforth was inspired by the 
beautiful work of Altman and Auerbach [10]. Let the 
operator Hit denote the creation operator of a hole pair on 
a plaquettc. The difference between pairs in the cuprates 
and in conventional BCS superconductors is that within 
the cuprates the pair size is very small which allows us to 
consider a reduction to plaquettes. This is also seemed 
to be supported by numerical calculations. 

A. Mapping To An XXZ Model 

Employing the Matsubara-Matsuda transformation 
[11] (trivially valid in any dimension) 

si = □!, 

Si = (18) 

we may map the most general Hamiltonian describing 
pair motions and interactions of all ranges 

Heff = -^2 tij ( a l D J + a Pi) 
(ij) 

+ J2v ij n inj + J' ^ □!□; + ..., (19) 

{ij) 1 
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FIG. 1. A plaquette state. Following Altman and Auer- 
bach's approach we tile the plane with non-overlapping pla- 
quettes. The small circles represent single holes. In the above 
we consider holes forming pairs on plaquettes. 

to a two dimensional XXZ model 

+ j zi3 s z s z - j' Y^ s t ) 2 

(ij) * 

-Y,h'Sl + ... (20) 

i 

with = 2t ij , J zij = V' 13 . 

Examining the quadratic (0(S 2 )) sector for nearest 
neighbors (\i — j\ = 1), we arrive at the standard nearest 
neighbor XXZ model. The XXZ model has two distinct 
phases: Within the Ising limit (large \ J Z \) the sponta- 
neous breaking of Z 2 symmetry implies phase separation 
of the bosonic hole pairs. By contrast, for large | J |, the 
system is in the XY limit and the ground state is a su- 
pcrfiuid. At the transition between these two phases, the 
model is SU (2) symmetric. Z 2 order (number order) is 
dual and conjugate to superconducting order (XY order). 
We have derived stripes in an earlier article [1] assuming 
weak staggered boundary conditions on the real S z (a Z 2 
symmetry order of yet another origin). 

B. Phase Separation or Not 

As to be expected, whether phase separation of hole 
pairs will transpire is determined by the sign of the 
charge-charge interactions Vij . 

If the charge-charge interaction is repulsive for all \i— j\ 
(i.e. Vij > 0), then in the extended spin model we will 
find exchange constants Jfj > 0. A positive J z (a repul- 
sive potential Vij) favors charge neutrality (S^ ot ) when- 
ever the background charge is taken into account. For an 
attractive potential Vij < leading to J Z j < 0, there is 
viable ferromagnetic order. In such an instance, the po- 



tential Vij allows (and favors) regions of different phase 
densities (i.e. phase separation). 

C. Superconducting Correlations, vortices, And A 
Numerical Value For Optimal Doping For Attractive 
Interactions 

Whenever the expectation value (5 2 ) vanishes we im- 
mediately find half a pair per plaquette or, equivalcntly, 
a doping of quarter of a hole per site. In order to make 
the XXZ model maximally superconducting, we would 
like to make the order as XY like as possible to avoid 
phase separation into bosonic hole pairs. If we indeed 
impose from the outset the absence of Z 2 like order for 
the number operator (S z ) we find a quarter of a hole per 
site 

^optimal = 1/4. (21) 

The largest S z (or number) fluctuations ((S z — (S z )) 2 } 
occur about the symmetric point (S z ) = 0. This fea- 
ture is far more generic than assuming a special model of 
two dimensional hard core bosons (as we have above). In 
high enough dimensions, the maximal fluctuations at the 
symmetry points may diverge signaling a critical point. 
Large number fluctuations enable well defined supercon- 
ducting phase (XY) order. A trivial bound on XY order 
reads 

(Si) < 5(5+1) -(S z ) 2 . (22) 

As the number operator is the dual disorder operator for 
the phase, in order to have maximal phase ordering, we 
wish to avoid vortex like number expectation values. A 
finite number density (n — 1/2), is analogous to a finite 
vortex density. This, in turn, is equivalent to an imposed 
external background magnetic flux density in a Joseph- 
son junction array. Non-uniformities in the charge order 
act as frustrating vortices. For maximal phase ordering 
we need to minimize the appearance of vortices. This 
entails lowering the effective magnetic flux background 
that spontaneously generated vortices must cancel. 

Another way of viewing matters at very low finite tem- 
perature (like a low temperature classical model) is as 
follows: If J z were zero then adding an additional field 
h! would corrupt the XY order- the previous classical 
ferromagnetic XY ground state- S± = (cos <f), sin <fi) with 
the phase <j) uniform over the entire plane. Now imagine 
turning on J z and making it negative (attractive interac- 
tions amongst pairs). This will only further corrupt the 
XY order. 

The relatively large deviation of the optimal doping 
found in the cuprates from 1 /4 suggests that we may not 
easily view these doped Mott insulators at optimal or 
near optimal doping as merely composed of pairs of the 
plaquette size with attractive interactions. 
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FIG. 2. The DMRG calculation of the t-J model by White 
and Scalapino. This and similar calculations typically employ 
open or cylindrical boundary conditions. 

V. PLAQUETTE STATES ON STRIPES 

The stripe constitutes a small system coupled to a large 
bath (the ambient antiferromagnet) from which it can 
draw holes. Depending on the chemical potential of the 
stripe, different hole densities will result. For ease we 
will invoke a simplifying yet non-essential assumption: by 
virtue of the large kinetic scale, the hopping amplitude t 
is much larger than all potential effects, i.e. Y^j Vij an d 
the magnetic alleviation energy (the energy gained by 
removing bad ferromagnetic bonds along the seam of the 
stripe) J. With such assumptions, the chemical potential 
trivially vanishes. In reality, due to the finite corrections 
of both effects we expect the chemical potential jj,/t ~ 0. 

As we will show assuming the constituents of the 
stripes to be of different geometries and nature (e.g. hole 
pairs on non overlapping plaquettes, diagonal or parallel 
hole pairs) different filling fractions will be found at the 
approximate particle-hole symmetric point fj, = 0. 

We will now require that every two consecutive rungs 
of the ladder form a unique plaquette (such that these 
plaquettes do not overlap) and consider all possible pair 
states within these plaquettes. 

As we will now show, if the chemical potential is indeed 
small by comparison to the kinetic hopping amplitude t, 
then assuming pair states in non-overlapping plaquettes 
will lead to the "correct" filling fraction. By "correct", 
we allude to the coincidence, within the large fi/t limit, 
with the filling fraction attained from an original large U 
Hubbard Hamiltonian for single fcrmions augmented by 
all possible hoppings and interactions amongst the holes. 
As we show in section(VII B), if we start from the original 
Hubbard model and allow holes to do anything that they 
wish from there (e.g. remain single, pair, phase separate) 
we will find that the net hole filling fraction is 1/4. This 
coincidence suggests that plaquette pair states are indeed 
viable candidates for the pairs states populating bond 
centered stripes. Moreover, both numerical [2], and mean 
field [12] calculations lead to quarter filling as well as 
experimental data [13]. 



These plaquette pair states are indeed the states ob- 
served numerically. Moreover, as emphasized by [10] they 
naturally display d- wave symmetry. 

As before, we may use the operator to denote the 
creation operator of a hole pair on a plaquette and sub- 
sequently employ the Matsubra-Matsuda transformation 
with no change to obtain once again Eqn.(29). 

As seen in Fig. (2), adding a hole pair along the bond 
centered stripe removes energetic ferromagnetic bonds. 
Such an effect will only serve as a weak sink favoring 
more hole pairs to drift into the stripe. 

If the charge-charge interaction, 

E VijUiUj - 2 J iij 

ij i 

= E km V2)K - 1/2) +5>E v * 

ij i j 

-2jJ2m, (23) 

i 

is repulsive for all \i — j\ (i.e. Vij > 0), then in the 
extended spin model we will find exchange constant J?- > 
for all \i — j\. A positive J z favors charge neutrality 
(5 t z ot ) when the background charge of (-1/2) (if ^ . Vij = 
2 J) is taken into account. For an attractive potential 
Vij < leading to J?- < 0, there is viable ferromagnetic 
order. The potential Vij allows (and favors) regions of 
different phase densities (i.e. phase separation). 

In the presence of a vanishing electronic chemical po- 
tential 

n = E v v - 2J > ( 24 ) 

3 

or, effectively, an external magnetic field h of a similar 
magnitude, there is no spontaneous symmetry breaking 
in the XXZ chain 

(S z ) = (25) 

and there is an average of one pair per two plaquettes 
or a quarter of a hole per site. This is indeed in accord 
with experimental observations [13] and the numerical 
calculations of White and Scalapino [2]. In reality, the 
chemical potential albeit small compared to t is not zero 
and a small deviation from 1/4 will be found. Nothing 
is inserted by hand here- the occupancy is dictated by 
the chemical potential (or magnetic field) which is finite 
yet may small compared to t by virtue of the various pa- 
rameters that it includes (i.e. J/t and j ^\ Vij). At its 
very core, the unbroken spin rotational symmetry cor- 
responds to particle-hole symmetry trivially present at 
fi/t = 0. This symmetry is lightly lifted by a small finite 
chemical potential shift relative to its vanishing value the 
particle- hole symmetric point. 

We will later on demonstrate that if we assume only 
symmetrized rung states for each of the individual holes 
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then no matter how the holes interact and whether they 
form pairs, phase separate, or not, a bond centered stripe 
must be nearly quarter filled when the chemical potential 
for the insertion of a hole is very small by comparison to 
the kinetic hopping scale t. 

To conclude, in the final analysis, we expect 

Sstripe w 1/4, (26) 

where the approximation sign signifies the relatively 
small (compared to t) chemical potential shift due to 
magnetic alleviation energies and hole-hole interactions. 

VI. NON PLAQUETTE PAIR STATES 

Here we show that unless the fundamental pair build- 
ing blocks are not chosen properly then the hole density 
may come out to be incorrect (when contrasted with the 
very general single hole problem for consistency in the 
large t limit). The filling fraction will be larger than a 
quarter found for the general hole problem before assum- 
ing anything about the possibility of pairing, phase sepa- 
ration etc. amongst holes. We will now analyze diagonal 
pairs which may be symmetrized and anti-symmetrized. 
Apart from some trivial numerical modifications, the con- 
siderations that we present here may be repeated word for 
word for "vertical" pair states extending along the rungs 
or "horizontal" pairs extending along the legs. There is 
a viable link between the correct symmetry of the pair 
state and its extension over different non-overlapping pla- 
quettes and the anticipated hole density. In order to ex- 
hibit this link and point to the viability of plaquctte pair 
physics, we will purposefully, focus on a pair state ansatz 
that will give an incorrect filling fraction in the \xjt = 
limit. More complicated and comprehensive treatments 
may employ much of what is known about the two leg lad- 
ders (for instance, the SO(8) symmetry present in their 
low coupling limit [14]). 

A. Diagonal Pairs 

We will now write down an effective one dimen- 
sional model for diagonal pairs of two possible chiralitics 
(left/right tilting) as shown in the right and left hand 
panels of Fig. (3). To give the reader a flavor of where we 
are heading the general argument may be summarized as 
follows- we will effectively make an exact projection of 
the two rung ladder along its axis to produce an effective 
one dimensional model. Henceforth we will denote the 
two natural diagonal pair chiralities (or polarizations) by 
a Greek index a = 1, 2 and further employ the shorthand 

A («=i) _ c t t 

— L 2,i,cr c l,i+l,-<7 

A 2 (a=1)t = ci, i+ i,_ <T c 2li , <r . (27) 



for the hole pair annihilation and creation operators of 
the diagonal polarization of Fig. (3) (a = 1) in terms of 
the electronic operators. In Eqn.(27) the first subscript 
(which is always one or two) marks the location of the 
hole on the ladder- on which of the two rungs it is found 
(upper or lower) while the second (i) denotes the location 
of the hole along the stripe (ladder) axis. The pair cre- 
ation operator removes two electrons and the pair annihi- 
lation operator creates two electrons. Note that in these 
operators, there is no spin index symmetrization/anti- 
symmetrization. The spin polarizations are dictated by 
the location along the ladder (whether or not i is even. 
We may conform to the convention that the spin sub- 
script a =t if i is even and that a —I along the odd 
rungs (i.e. we will choose our origin i = in such a way 
that the spin polarization of both electrons on that rung 
is a =t). The introduction of each diagonal pair states 
removes two bad magnetic bonds along the two rungs 
that it occupies. 

The number associated with these hole ( "h" ) pair op- 
erators is, as usual, 

»$ = [A^]tA^. (28) 

We therefore can write the most general effective one 
dimensional Hamiltonian 

^/ = -E*W )tA 5 w + A 5 wtA * (a) ) 

(ij) 

(ij) i,a 

where the ellipsis denote higher order terms in {A^}. All 
potential terms appear here including ones with infinite 
long range interactions (e.g. Coulomb effects) in which 
terms like decay as slowly as desired as a function of 
the separation \i — j\. Note that as a single hole cannot 
be doubly occupied (or an electron of fixed spin cannot 
be twice removed), no two diagonal bonds may share the 
same site; this may formulated in terms of an effective 
infinite hard core repulsion: V^"}^ =2 — ^2,'i +1 = 00 • 

Employing the Matsubara-Matsuda transformation 
once again we arrive at a flavored version of Eqn.(29) 
where we dress the various spins by the polarization (di- 
agonal orientation of the pair that they represent). Inso- 
far as commutators are concerned, the two polarization 
sectors are completely decoupled. Setting, for each a, 
A, = Sr and At = we see that the pair creation- 
and annihilation-operators form an SU(2) spin algebra. 
The effective Hamiltonian is a trivially flavored version 
of Eqn.(20): 

H ef f = -J2 J^iSf^sfV + Sf a) sf ] ) 

m 

(ij) i,a 
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FIG. 3. A hole pair advancing along the stripe. If we want 
to account for such motions (either as shown above, or result- 
ing from a higher order hopping process taking the pair in 
and out of the antiferromagnet surrounding the stripe back 
to the stripe), we must allow for overlapping plaquettes. As 
we will see, including such possibilities leads to higher filling 
fractions (i.e. 8 > 1/3 for a nearly vanishing relative chemical 
potential (fi/t)). 



-(«) 



(30) 



where the ellipsis denote higher order spin interactions. 
As before, to quadratic order, comparison with Eqn.(29) 
shows that Eqn.(30) is nothing but an XXZ model with 

and h = (2V + J) (with J the 



J 



±ij _ 



a/3 



2+y 



J af3 — 



V af3 



if V 



yn 



is constant for 



magnetic alleviation energy), 
all nearest neighbor links \i — j\ = 1. The only limita- 
tion of Eqs.(29, 30) is that they only allow for diagonal 
pairs as the fundamental building blocks (e.g. no single 
electron excitations are accounted for). 



B. Complete Fock Space (Lorentz) Algebra 

With an eye to things to come let reshuffle the two 
diagonal sectors a = 1 and a = 2 (corresponding to the 
left and right most real space configurations of Fig. (3). 

By a simple linear transformation 



S = S (a=1) + 5 (q=2) 



(31) 



the SU (2) ® SU (2) algebra of the two diagonal sectors 
transforms into a Lorentz (SO (A)) algebra. Physically, 



Ai 



A- +A- 
A! - A? 



(32) 



are the creation operators for the symmetric and anti- 
symmetric rung pair states. The symmetric combination, 
the angular momentum vector S (or L), forms an 5*0(3) 
subgroup: [5*, S J '] = i€ijkS k with the indices i,j and 
k now denoting the standard spatial components. The 
z-component 



5 : 



1 



E 



» 



(33) 



may assume the eigenvalues {1,0,-1} corresponding to 
a 2 x 2 patch of the ladder being electronically full, hav- 
ing one diagonal hole pair, and two having diagonal hole 

pairs. The symmetric combinations {Si} satisfy 5*0(3) 
algebra with spin one. The net number of electrons is 
given by the total magnetization of the symmetric (S* ot ) 
component. 



C. Restricted Symmetric Hole Pair State Basis 

We will now make the physical assumption (fortified 
by our earlier analysis [1]) that only the the low energy 
symmetric hole pair bonding states appear at low ener- 
gies. Defining, within the restricted symmetric hole pair 
state basis, 



57 



1 

75 



(A, 1 



A,; 



1 



b * " 2 



v 2 h 



[1 ^f^i-.syram] j 



(34) 



we see that 5 satisfies 5[/(2) algebra with 5=1/2 within 
the restricted symmetric bonding state space IIi(Aj ) Ti |0) 
with Ti = or 1, and where |0) denotes the vacuum state 
in which no holes are present. The bosonic hole pair 
states fill up the ladder as the electronic orbital states do 
the in atomic and molecular systems. We may regard this 
as an effective Hund's rule. At low hole filling fractions, 
only the low energy bonding states will be occupied. At 
high hole occupancies (in excess of half a hole per unit 
site), anti-symmetric anti-boding states will also appear 
and the full 50(4) algebra of the full blown unrestricted 
Fock space will raise its head. 

Substituting Eqn. (32,33,34) into Eqn.(29) and omit- 
ting any terms containing anti-bonding operators (which 
will take us out of the restricted symmetric pair state 
subspace) an arbitrarily high polynomial in {Si} will re- 
sult identical to Eqn. (20). The hard core repulsion term 
V[n 



(a=l) (a=2) 



(a=2) (a=l) 



transforms into an analo- 
gous hard core term VS*Sf +1 for the z-component of the 
spins (the distance between up spins cannot be smaller 
than two). 

Such a hard core term may be interpreted from first 
principles. A state containing nearest neighbor symmet- 
ric bonding pairs is not normalized to one. To have nor- 
malization we must consider the correlated four hole state 
2-1/2 E a =i,2 A" tA^lO). The reduced norm (probabil- 
ity) of the nearest neighbor bonding pairs vis a vis other 
bonding pairs is similar to having a large potential bar- 
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D. One Third Doping And Beyond 

Sans the hard core constraint, the model has its ground 
state at half filling (i.e. the density of up spins is a half): 
Whenever an XXZ model has a nearest neighbor hard 
core constraint (whenever the distance between up spins 
cannot be smaller than two), the density of up spins is 
equal to a third within the ground state [15]. A careful 
counting of the various diagonal pair states reveals that 
the hole occupancy within the stripe is expected to be 
1/3 in the large t limit. Repeating the same exercise for 
hole pairs along rungs and legs we arrive at similar large 
hole densities. 

VII. SINGLE ELECTRONIC DESCRIPTION 

A. Symmetric Rung States 

In an earlier work [1] , we have numerically obtained the 
single hole and hole pair energy spectrum on a stripe. Let 
us first restrict attention to the single electronic descrip- 
tion. The states along the various rungs form resonant 
bonding states. We will now make the physical assump- 
tion (fortified by numerics) that only the the low energy 
symmetric (or anti-symmetric) hole rung bonding states 
appear at low energies. We will now consider the sym- 
metric (anti-symmetric) smearing of a hole along the two 
legs of the rung. Defining, within the restricted symmet- 
ric hole rung state basis, the two operators 

Ci , ± ^^r l) ±c^ 2 \ (35) 

with c the electronic annihilation operator, a = 1, 2 now 
a leg index, and i a rung label, it is readily verified that 
these operators satisfy disjoint canonical Fermi anti com- 
mutation relations: {ci ; ±,c+ ± } = Sij. 

We will assume that the even parity symmetric bond- 
ing hole states are lower in energy and restrict ourselves 
to that basis. In the problem of physical relevance, 
\t^Zj'^~ 2 \ ^> all other hopping amplitudes and the hole 
will quickly resonate between the two rungs of the ladder 
before inching its way along the ladder axis. 

Once again, the hole states fill up the ladder as the 
electronic orbital states do the in atomic and molecular 
systems. We may regard this as an effective Hund rule. 
At low hole filling fractions, only the low energy bonding 
states will be occupied. At high hole occupancies (in ex- 
cess of half a hole per rung of a bond centered stripe or a 
hole per plaquette in the cuprate plane), anti-symmetric 
anti-boding states will also appear and the algebra of the 
full blown unrestricted Fock space will raise its head. 

In what will follow shortly, we will treat the symmet- 
ric rung holes as spinless fermions. An average density 



of half a spinless fcrmion (the particle-hole symmetric 
point) per rung corresponds to a quarter hole per site. 

We will consider the extended Hubbard model in the 
restricted basis of these low lying symmetric rung states. 

B. The Ideal Spinless Fermi Gas 

As we review in Appendix(A), within the general 
Hamiltonian Eqn.(l) in the limit of infinite on-site re- 
pulsion U, the charge degrees of freedom transform into 
those of a spinless Fermi system. In the bare Hubbard 
model, at U = oo the system reduces exactly an ideal 
(non-interacting) Fermi gas with a dispersion 

efe = — 2tcosk. (36) 

For a non-interacting Fermi gas, the zero temperature 
occupancy 

(n k ) =6( f i-e k ). (37) 

If inserting a hole leads to no change in the energy bal- 
ance then the chemical potential, by its very definition, 
vanishes. The energetics of the chemical potential acts 
as a Lagrange multiplier enforcing a certain average oc- 
cupancy. The energy of adding or removing an electron 
is the same (zero) at a chemical potential fi = 0. By 
particle-hole symmetry, at /i = the chemical potential 
lies in the middle of the band and the occupancy of the 
rung symmetrized spinless fermi particles is a half- and 
that of the holes is 1/4. The particle-hole symmetry of 
the kinetic energy band is dictated by the hermiticity of 
the Hamiltonian. In fact, even if the Hamiltonian in- 
cluded arbitrary long range hoppings {t r } and reduced 
to the general form 

H = -Y J t\i-iMcj + h.c), (38) 

ij 

in the infinite U limit, then by hermiticity whenever the 
chemical potential is zero the Hamiltonian will be particle 
hole symmetric. In that case, the density of rung sym- 
metrized spinless fermions will be a half, and the density 
of holes along the bond centered stripe will be 1/4. 

If the addition of holes (or removal of electrons) leads 
to a reduction in the magnetic strain energy along the 
stripe then the chemical potential for the electrons is re- 
duced (or, equivalently, that for holes, is increased). A 
removal of a bad bond along the rung leads to a lower- 
ing of the magnetic energy by J. Thus the hole chemical 
potential (i = J and the net, rung symmetrized, charge 
occupancy is given by an integral of (rife). The hole den- 
sity is ((rife/2)). This expectation value (hole occupancy) 
changes slightly with temperature according to the evo- 
lution of the Fermi function. 
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It is worthwhile viewing this, one last time, in terms of 
unbroken spin rotational symmetry in an effective XX(Z) 
model. 

We now apply the Jordan- Wigner transformation in 
order to arrive at XX(Z) model. As well known, the 
infinite U Hubbard Hamiltonian trivially transforms into 

H = —2t' s ^2 l {^i,xSi+\,x + Si,ySi+l t y). (39) 

i 

An infinitesimal chemical potential leads to an additional 
small magnetic field coupling to a z . The absence of spin 
symmetry breaking in this trivial example (no a z interac- 
tion appears in H) leads to the conclusion that the band 
is half filled (or that the hole density is 1/4). 
This is immediately seen taking note of 

St = 1/2 - c\ Cl . (40) 

A vanishing (<j z ) for the XY Hamiltonian of Eqn.(39) 
leads to (nj) = 1/2. Related back to our original prob- 
lem, this implies a quarter empty electronic bond cen- 
tered stripe. This value is seen here to be dictated by an 
unbroken particle-hole symmetry. 

If additional interactions next nearest neighbor pieces 
of the Coulomb and other interactions Vij n i n j) 
and the magnetic energy alleviation energy J are added 
to the Hamiltonian then the Hamiltonian for the spin- 
less symmetrized holes readily transforms to a full blown 
XXZ model with 

H = —2t ^ ^ {Sj X Si + l tX + SiySi + l,y) 

i 

+ Ef(!- S?/2)(l - S*/2) - J - £(1 - Sf/2) 

{ij) * 

= J2i-J±(S + S- + h.c.) + J z StS*\ 

{ij) 

+ £>S?. (41) 

i 

We first note that if the electrostatic energy V^+i is 
such that it equals the magnetic alleviation energy J then 
by the absence of spontaneous symmetry breaking in one 
dimension, (S z ) = and the bond centered stripe is quar- 
ter filled. If, hypothetically, S z would spontaneously de- 
velop a nonzero magnetization then this would imply that 
even if in an ensemble of stripes, the overall density of 
holes would be a quarter on average, different stripes in 
the ensemble would exhibit two different densities about 
that mean. 

An imbalance between V and J leads to nonzero h and 
to a finite value of (S z ) ^ 0- a deviation from quarter 
filling on every stripe. In the above we took into account 
both the Coulomb and magnetic alleviation effects in one 
go. 



In a very nice paper by Nayak and Wilczek [3] it was 
observed that 1/4 filling of a doped chain is indeed pre- 
dicted from the Bethe ansatz solution of Licb and Wu 
[16]. Some of the symmetries of the 1/4 filled point were 
noted. Here we emphasize that at low chemical potential, 
quarter filling not only coincides with symmetries, but is, 
in fact, a rigorous outcome of symmetry considerations. 
Here we also note that a finite magnetic alleviation en- 
ergy of the holes amounts to a shift in the chemical po- 
tential. The approximate calculations employed earlier 
are not mandatory for the determination of the hole oc- 
cupancy within the ground state; in order to account for 
magnetic effects, we may simply set the hole chemical 
potential n = J, and subsequently integrate the known 
ideal spinless fermi number density (n e ) up to that energy 
to obtain an exact result. 

APPENDIX A: A REVIEW OF TRIVIAL SPIN 
CHARGE SEPARATION IN ANY DIMENSION 
IN THE LIMIT OF LARGE ON SITE REPULSION 

Here we review trivial spin-charge separation within 
any Doped Mott insulator (of arbitrary dimensions) in 
the limit of large on site repulsion energy U — > oo. Let 
tpcr!,... cr N {xi, ...,xn) be an electronic or "hole" eigenstate 
of the general Hubbard Hamiltonian augmented by all 
possible higher order interactions (these may be sparked 
by electron phonon terms, Coulomb repulsions etc.) and 
all possible range hopping amplitudes t,t',t", ... etc. 

Recall that, trivially, ip ai ,...a N ( x ii ---iXn), 

• is anti-symmetric by virtue of its electronic con- 
stituents: 1p airT:j (Xi,Xj) = —IpcTjViiXjjXi). 

Note furthermore that exactly at the U — ► oo limit, 

• ijj(xi = Xj) = irrespective of the spin indices 
(7i, . .. , on • The reason for this nodal behavior is trivial- 
opposite spin occupation is forbidden by a divergent on 
site penalty U, and parallel spin occupation is strictly 
forbidden by the Pauli principle. Also, at finite tempera- 
ture the probability for having any two particles occupy 
the same state strictly vanishes. 

If we invoke this nodal, hard core, condition, then we 
may now remove the on site Hubbard repulsion from 
the original Hubbard Hamiltonian to obtain a general 
translationally invariant spin independent Hamiltonian 
H. Herein lies the crux of the trivial spin charge separa- 
tion at infinite U. 

The most general solution to the Schrodingcr equation 
H\ip) = E\ip) satisfying the last nodal condition with a 
spin independent Hamiltonian H (all terms apart from 
the on site interaction U are spin independent) is 

tp<7 U ...,a N (xi, ...,x N ) = x(oi, ...,a N )W(x!, ...,x N ), (Al) 

and linear superpositions of such degenerate solutions, 
where W is a solution to the Schrodinger equation. Here 
the trivial spin-charge separation is manifest. 
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Let us now impose the Fermionic statistics condition 



r spin r charge 



(A2) 



spin 



for all i ^ j. Note that [P< 

is spin independent. As P^ orfle = -[i^ n 
the charge permutation operator also commutes with H. 
This implies the function W must satisfy 



as the Hamiltonian 

1 _ _pij 

spin ' 



P 



charge 



W(xi, ...,Xi, ...,xj, ...,x N ) = 



a 



charge 



W(x 



1 , . . . , Xj , . . . , 



(A3) 



The identity [P^ harge f = 1 implies that a% arge = ±1. 

In general, the energy of a nodcless function IF is 
expected to be lower. However, the charge degrees of 
freedom cannot be completely symmetric: if the spinor 
x(ci, ctjv) i s fermionic (completely anti-symmetric) 
then it must vanish identically for N > 3 particles for 
the spin 5 = 1/2 electrons. In any dimension, the ef- 
fective Hamiltonian can now be written in terms of the 
charge only spinless degrees of freedom. If the Hamilto- 
nian is particle-hole symmetric (as expected by hcrmitic- 
ity of the most general arbitrary range kinetic terms) 
then whenever the chemical potential vanishes (whenever 
the insertion of a hole and a particle both cost zero en- 
ergy) , the average ensemble density of the correct spinless 
charge degrees of freedom must be a half. If, for instance, 
in any dimension, all low energy physics could be cap- 
tured in terms of spinless charge pair degrees of freedom 
W = F{{4>{xi,Xj)}) then in terms of these pairs, the 
ground state of the system would be half occupied. 

The charge degrees of freedom encapsulated in W can 
be symmetric with respect to, at most, single pairs (and 
anti-symmetric within these pairs) ; a higher symmetry is 
ruled out by the impossibility of a spinor anti-symmetric 
in three and more spin 1/2 indices (spinons cannot be 
fermionic). 

Let us now regress to one dimension. Here, hard core 
bosons cannot be distinguished from spinless fermions 
and different symmetry states of W are in fact degenerate 
at the U = oo point. This degeneracy is lifted as t/U 
becomes arbitrarily small but finite. As U is extremely 
large but finite = \ip a .= a .(xi = xj)\ < \i\) ai =- aj {x,i = 
Xj)\ -C 1 and local singlet correlations are generated. 

It has indeed been established by one dimensional 
Bethe ansatz [17] that at infinite U, the function W is 
a Slater determinant of single particle momentum eigen- 
statcs. 
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